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Abstract: 

We consider W-boson scattering at high-energy e + e~ colliders and study one-loop log- 
arithmic electroweak corrections within the Standard Model assuming a light Higgs boson. 
We present explicit analytical results for W + W~ — > W + W~. Using the equivalent vector- 
boson approximation, we have implemented these corrections into a Monte Carlo program 
for the process e + e~ — > v e v e W + W~ . The quality of the equivalent vector-boson approxi- 
mation and of the logarithmic high-energy approximation for the electroweak corrections 
is discussed in detail. The impact of the radiative effects is quantitatively analysed. The 
corrections are negative and their size, typically of the order of 10%, increases with energy 
reaching up to —20% and —50% at the ILC and CLIC, respectively. 
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1 Introduction 



One of the foremost open questions in particle physics concerns the mechanism of elec- 
troweak symmetry breaking. Depending on its realization in nature, the understanding 
of this subtle mechanism must be approached in different ways. If the Standard Model 
(SM) with a light Higgs boson is realized, the Higgs boson can be directly produced and 
its properties investigated. If the Higgs boson is heavy or absent, a complementary ap- 
proach must be pursuit. In this case, information can be extracted from the scattering 
of longitudinally polarized gauge bosons. In fact, at high energies, longitudinal vector 
bosons unveil their origin as Goldstone bosons and, by virtue of the equivalence theorem, 
reflect the dynamics of electroweak symmetry breaking. 

Accordingly, vector-boson scattering (VBS) looks very much different in these different 
scenarios of electroweak symmetry breaking. In the SM with a light Higgs boson the gauge 
sector remains weakly interacting and the cross section can be reliably predicted within 
perturbation theory. In the alternative scenario, where a light Higgs boson is absent, 
perturbative unitarity is violated, and the longitudinal gauge bosons must become strongly 
interacting at high energies thus allowing for non-perturbative restoration of unitarity [1]. 

In the past years, a variety of models have been proposed to parametrize the strongly- 
interacting electroweak gauge sector, and to recover unitarity (see for instance Ref. [2] and 
references therein). The common prediction is an enhanced production of longitudinal 
gauge bosons. However, the phenomenological consequences as well as the new particles 
the various models provide can be sensibly different. They can be classified into two main 
groups [3]. In the most optimistic case, one could expect many new resonances at future 
colliders. In the less favorable scenario, the mass of any new particle could be much bigger 
than the energy scale probed at the planned accelerators. In this case, the indirect effect 
of such particles would only consist in a slight increase of the VBS event rate at high 
energy compared to the predictions of the SM with a light Higgs boson. 

Several studies have been performed in order to estimate the possible reach of the 
future lepton and hadron colliders [2-8]. The answers strongly depend on energy and 
luminosity parameters. Also a good control of the SM background can prove essential, 
particularly for the less favorable case. In this respect, the scattering of W bosons within 
the SM with a light Higgs boson constitutes an irreducible background to any new-physics 
signal pointing towards a strongly-interacting VBS regime. 

In this paper, we consider VBS within the SM with a light Higgs boson at the planned 
e + e~ colliders. For the International Linear Collider (ILC) [9] we consider a centre-of-mass 
(CM) energy = 1 TeV and for the Compact Linear Collider (CLIC) [10] ^/s = 3TeV. 
More precisely, we focus on the production of W-boson pairs plus neutrinos in the re- 
action e + e~ — > z/ c z/ c W + W~. For the projected luminosity L = lab^ 1 , the experimental 
collaborations will collect thousands of events coming from VBS in the high-energy do- 
main, where a possible strongly-interacting regime of the weak gauge sector could manifest 
itself. 

To match the envisaged statistical precision, the SM predictions have to be computed 
beyond lowest order. Indeed, in the very same high-energy region of interest, the elec- 
troweak radiative effects are enhanced by electroweak Sudakov logarithms [11-15], i.e. 
double and single logarithms of the ratio of the scattering energy over the vector-boson 
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mass (recent progress in the evaluation of electroweak Sudakov logarithms is discussed in 
Ref. [16]). At 0(a), these corrections can reach several tens of per cent, as confirmed by 
various analyses performed for different processes at lepton and hadron colliders [17-26]. 
Hence, in the case at hand, they have to be taken into account in order to search for small 
deviations between data and SM predictions that might appear as a signal of strongly 
interacting electroweak symmetry breaking. 

So far, the process e + e~ — > u e u e W + W~ has been computed at Born level, and found 
promising for investigating electroweak symmetry breaking at high invariant masses of the 
produced W-boson pairs [7,27,28]. The aim of our work is to study the contributions of the 
electroweak Sudakov logarithms, which represent the dominant electroweak corrections at 
high energies. 

For the calculation of the electroweak corrections we use the equivalent vector-boson 
approximation following the approach of Ref. [29]. Within this approximation we only 
consider virtual 0(a) corrections to the WW-scattering subprocess. These corrections are 
calculated using the method of Refs. [12, 13] in the high-energy logarithmic approximation. 
As in Refs. [12, 13], the virtual photonic corrections are split into a symmetric-electroweak 
and a purely electromagnetic part, which originate from above and below the electroweak 
scale, respectively. The former part is in practice obtained by setting the photon mass 
equal to in the photonic virtual corrections, and is infrared finite. The infrared singu- 
larities are contained in the purely-electromagnetic part and are cancelled when including 
soft-photon bremsstrahlung. Both the symmetric-electroweak and purely-electromagnetic 
parts are included in our analytical results. However the latter is omitted in our numerical 
studies since it strictly depends on the experimental setup. 

The paper is organized as follows. In Sect. 2 we define the process and give the setup 
for the numerical evaluation. In Sect. 3 we describe the strategy of the calculation and 
discuss the quality of the used approximations. Numerical results are presented in Sect. 4, 
and Sect. 5 contains the summary. Explicit analytical results are listed in the appendices. 

2 Process definition and numerical setup 

We consider the production of a W-boson pair plus two neutrinos in electron-positron 
collisions: 

e+e" -> v e v e W + W~. (2.1) 

This process contains the VBS subprocess W + W~ — > W + W~ generically described by the 
first Feynman diagram in Fig. 1. The WW-scattering signal is not the only contribution 
to the final state in (2.1). The irreducible background, exemplified by the second and 
third Feynman graphs in Fig. 1, is indeed sizeable and must be properly suppressed in 
order to enhance the VBS signal-to-background ratio. A set of appropriate kinematical 
cuts to be imposed is summarized below. 

In our analyses, we use the input values [30] 

M w = 80.403 GeV, M z = 91.1876 GeV, 

M H = 120 GeV, m t = 174.2 GeV, (2.2) 
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Figure 1: Feynman diagrams contributing to the process e + e~ — > i> e u e W + W~ . The first 
graph on the left is a generic representation of the WW-scattering signal. The remaining 
two diagrams are examples of irreducible background. 

for vector-boson, Higgs-boson, and top-quark masses. All other fermions are taken to be 
massless. The sine s w and cosine c w of the weak mixing angle are fixed by 

2 2 M W 

C W = 1 — S W = ~Jy[2' (2-3) 

z 

Moreover, we adopt the so called G M -scheme, which effectively includes higher-order con- 
tributions associated with the running of the electromagnetic coupling and the leading 
universal two-loop m t -dependent corrections in the definition of a. Using 

= 1.16637 x 10" 5 GeV" 2 , (2.4) 

we have 

a = V2G m M^4/tt= 1/132.38.... (2.5) 

In the following sections, we present results for a CM energy s/s = 1 TeV, which can 
be reached at the ILC, and for a/s = 3 TeV as is planned for CLIC. In both cases we 
assume an integrated luminosity L = 1 ab _1 . 

Based on the study of Ref. [7], we have implemented a general set of cuts, proper 
for ILC and CLIC analyses, defined as follows. In the listed cuts, the numbers outside 
parentheses are used for yfs = 1 TeV, those within parentheses for ^/s = 3 TeV. 

• We require a W-boson transverse momentum Pt(W ± ) > 100(200) GeV, since the 
production of longitudinal vector bosons, i.e. the VBS signal, is enhanced for large 
scattering angles and high energies. Moreover, this cut removes events dominated 
by t-channel photon exchange in subprocesses. 

• We require |cos6 l (W ± )| < 0.8, where ^(W 1 * 1 ) is the angle of the produced W 1 * 1 boson 
with respect to the incoming positron in the laboratory frame, since the VBS sig- 
nal is characterized by central W-boson production. This cut also removes events 
dominated by t-channel photon exchange in subprocesses. 
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• We require \y(W + ) — y(W~)| < 2, where ^(W 1 * 1 ) is the rapidity of the produced W* 1 
boson defined as y = 0.51n[(P + P L )/(P — P L )] and E and P L are the W-boson 
energy and component of the momentum along the beam axis, respectively. This 
additional angular cut ensures the production of central W bosons also in the CM 
frame of vector-boson scattering. 

• We require a transverse momentum of the W-boson pair P T (WW) > 40(50) GeV. 
This cut suppresses the reducible background coming from the process e + e~ — > 
e + e~W + W~, i.e. 77 fusion, when the two produced e ± are emitted for- 
ward/backward. 

• We require a neutrino-pair invariant mass M(u c u c ) > 150(200) GeV. This cut re- 
moves W + W~Z production events, in which the neutrinos come from the Z-boson 
decay (see the last graph in Fig. 1). 

• We require a diboson invariant mass M(WW) > 400(700) GeV. Selecting high 
diboson CM energies allows one to test a possible strongly-interacting regime of the 
electro weak gauge sector. 

3 Strategy of the calculation 

In this section, we describe the main ingredients of our calculation. We summarize the 
adopted approximations and discuss their domain of applicability. 

We consider the process (2.1), with two on-shell W bosons and two neutrinos in the 
final state. For this process, exact lowest-order matrix elements are employed in our Monte 
Carlo, simultaneously accounting for signal and irreducible background. We moreover use 
the complete four-particle phase space and exact kinematics. 

Computing 0(a) electroweak corrections in leading-pole approximation, as in 
Refs. [31-33] and references therein, has revealed successful for analysing WW physics 
at LEP2. A similar philosophy can be adopted for the incoming bosons in VBS at ener- 
gies that are large compared to the gauge-boson masses, since this process is dominated 
by small invariant masses of these incoming bosons, which are thus relatively close to 
their mass shell, even though these invariant masses are actually negative. 

We thus compute the 0(a) electroweak corrections to the process (2.1) in equivalent- 
vector-boson approximation (EVBA). As discussed in the introduction, we do not include 
real photonic corrections. For the virtual corrections we work in the logarithmic approx- 
imation and we restrict our calculation to the infrared-finite part coming from above the 
electroweak scale. These corrections correspond to the case where the photon has ef- 
fectively the mass Myy and are precisely defined in Ref. [12]. This approach is sensible 
since at high energies, the electroweak corrections are dominated by double and single 
logarithms of the ratio of the energy to the electroweak scale. Hence, keeping only the 
terms proportional to aln 2 (s/M{y) and aln(s/M^), where s is the CM energy of the 
VBS subprocess, provides the bulk of the radiative corrections. 
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3.1 Equivalent vector-boson approximation 

In EVBA, the process e + e~ — > v c v c W + W~ is entirely described by the subset of Feynman 
diagrams generically represented by the first graph in Fig. 1. The approximation considers 
in fact only those contributions to the final state that come from the scattering of the two 
vector bosons emitted by the incoming particles. The goodness of the approximation thus 
relies on the assumption that such contributions are indeed the dominating ones in the 
considered kinematical domain. This hypothesis depends on the process at hand, and can 
only be checked against an exact computation. In this section, we discuss the reliability 
of the EVBA in describing the process e + e~ — > z/ e z/ e W + W" and its validity domain. 

In implementing the EVBA, we follow the approach of Ref. [29]. This method pre- 
serves the exact kinematics of the process, a very useful property for imposing realistic 
cuts. 

We assign the following set of momenta Pi and helicities Aj = 0, ±1 to the particles 
involved in the process we are considering: 

e + (pi, +) e~(p 2 , -) -> jy c (ps, -) v c (p 4 , +) W + (p 5 , A 5 ) W~(p 6 , A 6 ). (3.1) 

In EVBA only left-handed lepton chiralities are relevant, since the electrons and neutrinos 
couple always to W bosons. 1 Using the unitary gauge and writing the propagators of the 
two incoming W bosons, emitted by the initial e ± and exchanged in t channel (see Fig. 1), 
as a sum over the vector-boson polarizations 



f-Ml V Ml) X= % A pi-Mi 

the exact amplitude corresponding to the first graph in Fig. 1 assumes the form 

1 1 



^e+e-^ePeW+W" {pi ? ^ pg . 



x £ M c+ ^ w+ ( Pl ,p 4 ,q + ; A + ) .M e -^ w -(p 2 ,p 3 ,<Z-;A_) 

A+,A_=-1,0,1 

x ^w + w^w + w- (g+) q _,p 5 ,p 6] A + , A_, A 5 , A 6 ). (3.3) 

In EVBA the off-shell amplitude for WW scattering, A4 W+W ^ w+w ; j s replaced by 
a suitably defined on-shell amplitude. In this way gauge invariance of this amplitude is 
ensured, and artifacts from using an incomplete off-shell amplitude are avoided. Modi- 
fications of the on-shell amplitude are necessary in order to describe the dependence of 
the off-shell amplitude on the off-shell masses q\ to a satisfactory accuracy. We here 
follow the approach of Ref. [29] which describes the extrapolation to off-shell masses by 
simple proportionality factors for each incoming vector boson. These are chosen to be 
equal to 1 for transverse bosons. For each longitudinal W + or W~ boson, the on-shell 



1 Instead, the exact matrix elements that we employ for our tree-level predictions receive (background) 
contributions also from right-handed leptons. 
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OBorn(e + e -> z/ c z/ e W + W ) 


[TeV] 


O-cxact [fb] 


0"EVBA [fb] 


Aevba [%] 


1 


0.595 


0.479 


19.5 


3 


3.507 


3.471 


1.0 



Table 1: Exact lowest-order cross section (second column) as well as total cross section in 
EVBA (third column) and their difference in per cent of the exact result (fourth column). 
Kinematical cuts as in Sect. 2 are applied. 



amplitude is multiplied by a factor / ' \J —q± in order to describe the singular behaviour 
e \=o(Q±) ~ 1/ \J — Q± °f the off-shell longitudinal polarization vectors. One can thus write 

1 1 



M e e 



EVBA 



X 



+ y P c W+W- 



(Pl,P2,P3,P4,P5,P6; A 5 , A 6 ) 



w 



c+^P W+ 



(Pi,P4,9+; A+) M c 



(P2,P3,5-; A_) 



X ^W + W-^W + W- (g o„ ; g 0„^ 5;p6 . A+? Ag) 



X 



Mi 



w 



<V,0 + 5 A+,± 



M 



w 



(3.4) 



where q°± are the on-shell projected momenta of the two incoming W bosons. The defini- 
tion of the on-shell projection is given in Appendix A. Note that the W-boson momenta 
in the matrix elements M e+ ^ VcSN+ and Ai° ^ cW are not projected on shell. 

In order to proceed, in Ref. [29] the squared amplitude was considered and the contri- 
butions of the matrix elements yvi e+ ^^w + anc ^ j^e ^ c w were transformed into vector- 
boson luminosities. Instead, we work at the matrix element level and compute the three 
amplitudes on the right-hand side of (3.4) with the help of PHACT [34], a routine based on 
the helicity-amplitude method of Ref. [35]. 

In Table 1, we show the comparison between the EVBA and the exact lowest-order 
result for the total cross section. We select the kinematical domain where we expect 
new-physics effects related to strongly interacting vector bosons to be enhanced. Such a 
region, characterized by high diboson invariant masses and large scattering angles of the 
two produced W bosons, is selected via appropriate cuts described in detail in Sect. 2. 
At y/s = 1 TeV and \fs = 3 TeV the accuracy of the EVBA for the total cross section, 
Aevba = (o"cxact — 0"EVBA)/ocxact, amounts to about 20% and 1%, respectively. 

The integrated cross section gives only a partial information on the goodness of the 
EVBA. In order to display more extensively the reliability of this approximation in the 
selected kinematical domain, in Figs. 2 and 3 we analyse distributions in both energy-like 
and angular-like variables. In particular, we consider four observables of interest: 



diboson invariant mass M(WW), 
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dAf(WW) [P b /GcV] 



r(WW) [P b / GcV l 




Af(WW) [GeV] 




P T (WW) [GeV] 




[PVGcV] 




y(WW 



P T (W) [GeV] 



Figure 2: Lowest-order distributions for = 1 TeV from the exact matrix elements and 
in EVBA: invariant mass of the diboson pair (upper left), transverse momentum of the 
diboson pair (upper right), rapidity of the diboson pair (lower left), transverse momentum 
of the produced W boson (lower right). The inset plots show the relative difference Aevba 
in per cent. Standard cuts are applied. 



diboson transverse momentum Pr(WW), 
diboson rapidity y(WW) = 0.5 In 
W-boson transverse momentum Pp(W) = Pp(W + ). 



g(WW)+P L (WW) 
E(WW)-P L (WW) 



The distributions in Pr(W~) and Pr(W + ) are identical. We plot the lowest-order results 
of the EVBA and of the exact calculation for the two collider energies \fs = 1 TeV and 
y/s = 3 TeV in Fig. 2 and Fig. 3, respectively. The inset plots show the difference in 
per cent between the two results, i.e. A EV ba = (do" ex act — dcrEVBA)/dcr exact . For ILC and 
CLIC, in most of the regions that are not statistically irrelevant, the difference is below 
25% and 20%, respectively. The increase for small M(WW) is not problematic, since the 
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-(/(WW) P T (W) [GeV] 



Figure 3: Lowest-order distributions for = 3 TeV from the exact matrix elements and 
in EVBA. Same conventions as in Fig. 2. 

radiative corrections are small in this region (see Sect. 4). As a remark, let us add that 
the agreement between EVBA and exact result reached in this analysis should not be 
taken as the best one can do. The quality of the approximation is sensibly cut- dependent. 
In this paper, following Ref. [7] we have adopted generic cuts and, in principle, the EVBA 
could be improved imposing more stringent constraints. The choice of the set-up depends, 
however, on the particular analysis to be performed. 

In the following, we use the EVBA only for computing the 0(a) electroweak correc- 
tions. The lowest-order results are calculated exactly. The 0(a) inaccuracy associated 
with the EVBA can be estimated to be of the order of the product of the 0(a) correc- 
tions (computed in EVBA) times the inaccuracy of the EVBA at tree level. Since the 
radiative effects typically amount to 10-30%, the inaccuracy of the EVBA translates into 
a few-percent error for the full result. As shown in Sect. 4, apart for the case of very high 
P T (WW), the 0(a) uncertainty associated with the EVBA is always smaller than the 
expected statistical error at the ILC and CLIC. 
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e 



e+ 




Figure 4: Structure of the virtual factorizable corrections in EVBA with one-loop contri- 
butions in the blobs. 

3.2 Virtual electroweak corrections 

In analogy with the leading-pole approximation for virtual corrections, the EVBA allows 
two types of radiative contributions, factorizable and non-factorizable ones. The former 
are those that can be associated to either the emission of one of the two incoming W bosons 
from the beam particles or the VBS subprocess. The latter are those connecting these 
subprocesses. 

The non-factorizable corrections in leading-pole approximation consist only of pho- 
tonic contributions. These corrections have been evaluated for W-boson pair production 
in e + e~ annihilation in Refs. [36,37]. There it was found that all infrared and mass- 
singular logarithms cancel between virtual and real corrections and that the remaining 
effects are small. In EVBA non-factorizable corrections do not only result from photon 
exchange between the subprocesses of W-boson-pro duct ion and WW scattering but also 
from analogous exchanges of massive gauge bosons. These contributions deserve further 
investigations. In this paper, we do not consider the non-factorizable corrections but 
restrict ourselves to the calculation of the factorizable contributions. 

The virtual factorizable corrections are represented by the schematic diagram of 
Fig. 4, in which the big blobs contain all one-loop corrections to the incoming W-boson- 
production and on-shell WW-scattering subprocesses. The corresponding matrix element 
can be written as 



5M 



c+c-^ c P c W+ W" 
virt,EVBA,fact 




1 



1 



M 



c -^u c W x 



Born 



M 
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e+^P c W+ e-^i/ c W: W+ W7 W7 

+ - /W Born °- /W virt -^'Bori^on 

e+^P c W+ c--*u c W~ W+ W~ W7 

_|_ KA + KA fiM + 5 6 

^ JVl Bom JVl Born 0JVl virt,on 

M w 



X 



Mw X i r 

"9+ 



(3.5) 



where 5M. Virt + , <5.M virt , and SAi virt + on 5 6 denote the virtual correc- 

tions to the matrix elements for the W-boson-production and on-shell WW-scattering 
subprocesses. The index 'on' indicates that on-shell vector-boson momenta are used to 
calculate these matrix elements. 

We calculate the factorizable 0(a) virtual corrections in logarithmic high-energy ap- 
proximation including single and double enhanced logarithms, i.e. contributions propor- 
tional to a\n 2 (s/M^) and a ln(s/M^), where s is the CM energy of the scattering sub- 
process. The logarithmic approximation yields the dominant corrections as long as CM 
energies and scattering angles are large. Pure angular-dependent logarithms of the form 
am 2 (s/f) and aln(s/f), with r equal to the Mandelstam variables i and u of the WW- 
scattering subprocess, are not included. However, angular-dependent terms of the form 
aln(s/f) ln(s/M(^) are taken into account. The validity of the results relies therefore on 
the assumption that all invariants are large compared with and approximately of the 
same size 

s ~ |i| ~ |u| > M%f- (3.6) 

This implies that the produced gauge bosons should be energetic and emitted at suffi- 
ciently large angles with respect to the beam. This is precisely the kinematical region 
where effects due to a possible strongly interacting regime of the gauge sector are maxi- 
mally enhanced. In this region, the accuracy of the logarithmic high-energy approximation 
is expected to be of the order of a few per cent. We can thus reasonably adopt this ap- 
proximation at e + e~ colliders with energy in the 1-3 TeV range, where the experimental 
error is at the few-per-cent level. Since the emission subprocesses of the two incoming 
W bosons involve no large energy variable (they peak at \q±\ ~ -^w), the corresponding 
virtual corrections vanish in the logarithmic approximation. As a consequence, we do 
not consider the first two contributions on the right-hand side of (3.5) in the following. 
Moreover, for the W + W~ — > W + W~ subprocess we take into account only the corrections 
to the dominating channels involving four transverse (TTTT) or two transverse and two 
longitudinal (LLTT, LTLT, TLTL) gauge bosons. The contributions of the channels with 
an odd number of longitudinally polarized W bosons are suppressed by Myf/y/1, and 
those of the channels LTTL and TLLT by M^/s. Moreover, the configurations with two 
final-state longitudinal W bosons are numerically small within the SM with a light Higgs 
boson. As shown in Table 2, for Mh = 120 GeV, the cross section for the production of 
two longitudinal W bosons is suppressed by a factor 50 or 100 compared to the full result 
for y/s = 1 TeV or 3 TeV, respectively. 

The analytical expressions for the 0(a) virtual corrections to WW scattering in the 
high-energy limit are given in Appendix B. The formulas are rather compact and easy to 
implement. In our default set-up, we have used the version with exact SU(2)-transformed 



10 



OBom(e + e -> z/ e z/ e W|W A ,) 


[TeV] 


a TT [fb] 


a TL [fb] 


a LT [fb] 


a LL [fb] 


a tot [fb] 


1 


0.500 


0.0410 


0.0410 


0.0134 


0.595 


3 


3.111 


0.1786 


0.1786 


0.0390 


3.507 



Table 2: Born cross section for the process e + e _ — > z/ c z/ e W^W^, and various polarizations 
(A, A' = T, L) of the produced W bosons. Kinematical cuts as specified in Sect. 2 are 
applied. 

lowest-order matrix elements. This means that we use the complete expression as in 
(B.21). More precisely, we use (B.27), (B.30), and (B.33). We have verified that the 
numerical results based on (B.28), (B.31), and (B.34), which are obtained by making 
use of the high-energy approximation for the SU(2)-transformed Born amplitudes given 
in Appendix D, are in very good agreement. For all results shown in the following, the 
difference between the two methods is in fact at the per-mille level. This comparison 
confirms the reliability of the high-energy approximation for the Born matrix elements, 
under which the correction factor can be factorized and expressed in a very simple form, 
leading to considerable decrease in CPU time. 

In order to estimate the accuracy of the logarithmic high-energy approximation in com- 
puting the radiative effects, we have moreover compared our results with the complete 
0(a) corrections to the on-shell W A+ W A — > W^ 5 W^ 6 process [38]. For this comparison, 
we have included the electromagnetic terms given in Appendix E, evaluated for a soft- 
photon cutoff AE = 0.5y/§. In Fig. 5, we show the difference between the complete 0(a) 
cross section and the O{o) result in logarithmic high-energy approximation, normalized 
to the Born cross section, as a function of the CM energy. The three curves represent 
the three polarized cross-section differences for A + A_A 5 A 6 = TTTT, LLTT, TLTL, respec- 
tively. They all display a plateau at high energy, showing that the difference between 
approximate and exact result is due to terms that do not grow with energy in that region. 
This implies that the high-energy behaviour of the radiative effects is well reproduced by 
the adopted approximation which takes into account only energy-enhanced logarithmic 
terms. 

The terms, omitted in the logarithmic approximation might be, however, not negligi- 
ble. Figure 5 shows indeed that the agreement between exact and approximate 0(a) cross 
section is within a few per cent for the first two polarization configurations, while it goes 
up to about 13% for the third one. This discrepancy is unexpectedly large. It could be un- 
derstood by performing a complete high-energy approximation including non- logarithmic 
terms. However, for the unpolarized cross section this discrepancy is still tolerable. The 
lowest-order cross section is in fact dominated by the TTTT configuration, and receives 
only a 10% contribution from TLTL, as shown in Table 2. We can thus safely assume our 
0(a) inaccuracy to be of the order of a few per cent. 
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V 7 ! [GeV] 

Figure 5: WW-scattering subprocess W|" + W^"_ — > W^ g W^" 6 . Difference between complete 
0(a) cross section (a^ 111 ) and 0(a) cross section in logarithmic high-energy approximation 
(<r HE ), normalized to the lowest-order result. In the legend, T and L refer to the transverse 
and longitudinal polarizations of the two incoming and two outgoing W bosons from left 
to right. 

4 Numerical results 

In this section, we illustrate the effect of the logarithmic electroweak corrections on the 
production of a W + W~ pair plus missing energy at future e + e~ colliders. We consider 
the process (2.1) in the numerical setup given in Sect. 2. We analyse the behaviour of the 
VBS included in (2.1) over the kinematical region characterized by large diboson invariant 
masses and large scattering angles of the produced W bosons. This is the domain where 
possible new-physics effects entering the VBS would be maximally enhanced. 

We focus on the SM predictions for a light Higgs boson. As pointed out by Bagger 
et al. [4], in this case one would observe the production of mostly transversally polarized 
W bosons in the high WW invariant- mass region. The longitudinal spin configuration is 
in fact strongly suppressed in the presence of a light Higgs boson. For the specific case at 
hand, this is shown in Table 2 for Mh = 120 GeV and two possible setups. 

In Table 3, we show the impact of the corrections on the total cross section for 
e + e~ — > h> e u e W + W~ and compare it with the expected statistical error. The second 
column contains the lowest-order result. The third and fourth entries respectively display 
the 0(a)-corrected cross section and the contribution of the one-loop corrections relative 
to the Born result, Aew = (c — c r Born)/c r Bom- The fifth column provides an estimate of the 
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<r(e + e - 


-> i/ e ^W+W-; 






Vs [TeV] 


0"Born [iDj 


a [fb] 


A EW [%] 


A [%] 


l/V^Born [%] 


1 


0.595 


0.556 


-6.7 


1.3 


4.1 


3 


3.507 


2.897 


-17.4 


0.2 


1.7 



Table 3: Total lowest-order cross section (second column) as well as total 0(a) cross 
section (third column) and electroweak corrections in per cent of the lowest-order result 
(fourth column), including their uncertainty (fifth column). The last entry shows the 
statistical error for an integrated luminosity L = lab -1 . Kinematical cuts as in Sect. 2 
are applied. 

uncertainty of the one-loop contributions due to the EVBA, A = A E w x A E vba, which is 
obtained combining the one-loop corrections with the uncertainty of the EVBA at Born 
level. For the considered process, the electroweak radiative effects are negative and of the 
order of —5% to —20%. This has to be compared with the last column of Table 3, where 
we show an estimate of the statistical error based on an integrated luminosity L — 1 ab _1 . 
As can be seen, the electroweak corrections are quite important. Already comparable 
with the statistical uncertainty at the ILC, they further increase at higher energies giving 
rise to a 10a effect at CLIC. 

The influence of the 0(a) corrections is highly dependent on the cuts imposed and 
the selected kinematical domain. Also, distributions can be differently affected by the ra- 
diative corrections. We illustrate this point for the sample variables defined and analysed 
in Sect. 3.1. The 0(a) effects on these four observables at \fs = 1 TeV and 3 TeV are 
displayed in Fig. 6 and Fig. 7, respectively. The upper curves represent the lowest-order 
differential cross section, the lower ones the corresponding corrected result. 

We first consider the case of yfs = 1 TeV. In the left-upper plot of Fig. 6, we show 
the distribution in the diboson invariant mass. This variable, representing the CM energy 
of the VBS subprocess, gives direct access to the energy scale at which new physics 
could appear. In absence of a light Higgs boson, the SM VBS amplitudes would violate 
perturbative unitarity at high CM energies. In order to recover it, new physics should 
manifest itself at those scales. Hence, at future colliders it will be useful to analyse the 
diboson production (plus missing energy) at the highest possible M(WW) values. In this 
region, the O(a) corrections are enhanced. They can go up to —18%, as shown by the 
corresponding inset plot. The increase of the radiative effects with M(WW) is a typical 
effect of large logarithms of Sudakov type. The electroweak corrections should therefore 
be included to match the experimental accuracy. 

A second variable of interest, P T (WW), combines energy and angle information. This 
observable is expected to be sensitive to a strongly interacting VBS signal at low and 
intermediate values, say for Px(WW) between 50 and 300 GeV [7]. In this region, the 
radiative effects are of order —5% to —7%, thus smaller than in the case discussed above. 
Their size is, however, comparable with the statistical accuracy in the considered range. 
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Figure 6: Distributions for y^s = 1 TeV in lowest order and including logarithmic cor- 
rections: invariant mass of the diboson pair (upper left), transverse momentum of the 
diboson pair (upper right), rapidity of the diboson pair (lower left), transverse momen- 
tum of the produced W boson (lower right). The inset plots show the relative O{o) 
corrections A EW in per cent normalized to the lowest-order results. Standard cuts are 
applied. 

The same conclusion holds for the pure angular-like variable, 3/ (WW), shown in the left- 
lower plot, which gets corrections in the range between —7% and —10%. The transverse 
momentum Pt(W) on the right-lower plot displays instead a behaviour analogous to 
M(WW). Note that for all considered distributions the corrections reduce the tree- level 
SM predictions. 

In Fig. 7, we show the same set of distributions as above for ^/s = 3 TeV and the 
corresponding cuts specified in Sect. 2. As expected the impact of the radiative corrections 
increases with the collider energy since this allows for higher CM energies, which translate 
into higher diboson invariant masses and transverse momenta. This behaviour is well 
depicted in all four plots. For the distributions in Pr(WW) and y(WW) the corrections 
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Figure 7: Distributions for 3 TeV in lowest order and including logarithmic corrections. 
Same conventions as in Fig. 6. 

are in the range between —10% and —20%, apart from the region of very high Pt(WW) 
where statistics is small. For the M(WW) distribution the electroweak corrections grow 
from —10% to —50% with increasing invariant mass and for the Pt(W) distribution the 
effect is similar. For the distribution in y(WW), the corrections are large at the maximum 
of the distribution at small y(WW). For the other distributions, the cross sections get 
small where the corrections get large. Still, the 0(a) corrections are statistically relevant. 
This is illustrated in Table 4, where the size of the electroweak corrections Aew and the 
0(a) uncertainty A = Aew x Aevba, which results from the EVBA, are compared with 
the statistical accuracy A stat based on the projected luminosity L = 1 ab -1 . Dividing the 
range of energy-like variables in 200 GeV intervals, we find that the influence of the 0(a) 
corrections is not washed out by the binning. In the more statistically relevant bins, the 
radiative corrections ranging from —10% to —30% give rise to a 3-6<r effect. In the tails of 
the distributions, the 0(a) contributions can increase up to —50%, still being bigger than 
the estimated experimental accuracy. We also observe that, apart from the region of very 



15 



bin 


M(WW) 


bin 


P T (WW) 


[GeV] 


Ncvt 


Astat[%] 


A EW [%] 


A[%] 


[GeV] 


N cvt 


A sta t[%] 


A EW [%] 


A[%] 


700-900 


1719 


2.4 


-10.4 


0.7 


50-250 


1467 


2.6 


-15.9 


0.7 


900-1100 


916 


3.3 


-18.2 


1.2 


250-450 


1119 


3.0 


-16.6 


2.0 


1100-1300 


446 


4.7 


-25.1 


1.8 


450-650 


550 


4.3 


-17.0 


0.5 


1300-1500 


216 


6.8 


-30.4 


0.7 


650-850 


249 


6.3 


-21.4 


6.5 


1500-1700 


106 


9.7 


-36.8 


0.8 


850-1050 


92 


10.4 


-32.6 


24.9 


1700-1900 


53 


13.7 


-43.3 


1.5 












1900-2100 


26 


19.4 


-49.5 


2.2 












bin 


P T (W) 


bin 


y(WW) 


[GeV] 


iV evt 


A s tat [%] 


Aew[%] 


A[%] 




iV evt 


A stat [%] 


A EW [%] 


A[%] 


200-400 


1765 


2.4 


-13.4 


0.9 


0-0.25 


549 


4.3 


-21.1 


3.2 


400-600 


1140 


3.0 


-19.4 


1.2 


0.25-0.50 


486 


4.5 


-18.9 


1.1 


600-800 


424 


4.9 


-24.3 


3.6 


0.50-0.75 


374 


5.2 


-15.6 


1.4 


800-1000 


134 


8.6 


-27.8 


5.5 


0.75-1.00 


235 


6.5 


-12.4 


3.0 


1000-1200 


35 


16.8 


-28.6 


8.1 


1.00-1.25 


101 


9.9 


-9.9 


3.3 



Table 4: Impact of electroweak corrections on binned distributions at y/s = 3 TeV. The 
first and sixth columns show the bin. For each variable, the four entries from left to right 
give the total number of events for a luminosity L = 1 ab" 1 , the corresponding statistical 
accuracy, the size of the O(a) electroweak corrections relative to the Born result, and the 
estimated one-loop uncertainty due to the EVBA. 



high Pt(WW), the uncertainty associated with the EVBA never exceeds the statistical 
error. 

In Table 5 we give analogous results for y^s = 1 TeV. In this case, the radiative effects 
are comparable with the statistical accuracy. Their significance is strictly dependent on 
the binning. 

In our analysis, we have not included any option on the polarization of the initial 
beams. In the process e + e~ — > z/ c z/ e W + W~, the VBS signal is purely given by the W- 
boson scattering, as illustrated in the first Feynman diagram of Fig. 1. Only the right- 
left combination e^e^ thus contributes to the signal. The irreducible background (see 
for instance the last diagram in Fig. 1) can receive instead contributions also from other 
helicity configurations. The possibility of selecting a given initial polarization has thus two 
advantages. In first place, it helps in suppressing the background. As a second benefit, it 
increases the statistics. Assuming a polarization efficiency of 80% and 60% for electron 
and positron, respectively, we would have in fact an increase of about a factor two in the 
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bin 


M(WW) 


bin 


P T (WW) 


[GeV] 


N cvt 


A stat [%] 


A EW [%] 


A[%] 


[GeV] 


Ncvt 


A stat [%] 


A EW [%] 


A[%] 


400-500 


363 


5.2 


-3.5 


1.0 


40-140 


316 


5.6 


-6.1 


1.3 


500-600 


157 


8.0 


-9.4 


1.1 


140-240 


215 


6.8 


-7.4 


1.6 


600-700 


57 


13.2 


-16.0 


0.5 


240-340 


59 


13.0 


-7.6 


0.5 


700-800 


16 


24.8 


-19.9 


2.0 












bin 


P T (W) 


bin 


j/ (WW) 


[GeV] 


A^evt 


A sta t[%] 


Aew[%] 


A[%] 




N evt 


A sta t[%] 


A EW [%] 


A[%] 


100-200 


259 


6.2 


-4.8 


0.7 


0-0.25 


180 


7.4 


-6.5 


0.7 


200-300 


259 


6.2 


-7.9 


1.6 


0.25-0.50 


94 


10.3 


-6.9 


2.1 


300-400 


72 


11.8 


-9.4 


3.1 


0.50-0.75 


23 


20.8 


-8.3 


3.8 



Table 5: Impact of electroweak corrections on binned distributions at ^/s = 1 TeV. Same 
conventions as in Table 4. 



number of events at fixed luminosity. In this set-up, the relevance of the radiative effects 
would be further enhanced. 

5 Conclusions 

If the Higgs boson should be heavy or absent, the scattering of longitudinal electroweak 
gauge bosons can provide information on the mechanism of electroweak symmetry break- 
ing. An irreducible background to signals of new physics in vector-boson scattering is 
provided by the Standard Model contribution for a light Higgs boson. In order to be 
able to disentangle possible small new-physics effects a precise knowledge of the latter is 
required. 

We have studied electroweak radiative corrections to W + W~ scattering at high-energy 
e + e~ colliders. We have used the equivalent vector-boson approximation and calculated 
the factorizable one-loop corrections in the high-energy logarithmic approximation. Cor- 
rections to the splitting of the W bosons from the incoming particles as well as non- 
factorizable corrections have not been taken into account. We have presented explicit 
analytical results for the logarithmic corrections to W + W~ — > W + W~. These have been 
implemented into a Monte Carlo program for the process e + e~ — > u e i' e W + W~ together 
with the complete lowest-order matrix elements. 

We have defined a set of cuts suitable for the analysis of the scattering of strongly 
interacting W bosons and investigated our approximations within this setup. In kine- 
matical regions that are statistically relevant and receive non-negligible corrections, we 
find that the effective vector-boson approximation agrees with the complete lowest-order 
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prediction within about 20% to 25%. The logarithmic approximation reproduces the com- 
plete one-loop corrections for the unpolarized cross section of W + W" — > W + W" at the 
level of a few per cent. These approximations cause an uncertainty of our predictions at 
the level of a few per cent. This uncertainty is always comparable to the statistical error. 
In principle it can be further reduced by introducing appropriate cuts that improve the 
accuracy of the EVBA. But, in general, for a given set of cuts, a better precision can be 
achieved only by means of an exact calculation. 

The size of the electroweak corrections depends strongly on the cuts and the considered 
observable. We have studied their effect on the total cross section and four physically 
interesting distributions. Within our set-up, the corrections to the total cross section 
amount to —7% and —17% for y/s = 1 TeV and y/s = 3TeV, respectively. For the 
distributions in the transverse momentum and the rapidity of the W-boson pair they are 
of similar size. In the distributions in the invariant mass of the W-boson pair and the 
transverse momentum of a W boson, the corrections are negative, of the order of 10% and 
increase in magnitude with increasing energy. They can reach up to —20% and —50% 
for y/s = 1 TeV and y/s = 3 TeV, respectively. In summary, the electroweak corrections 
reduce the Standard Model predictions by a sizeable amount that is comparable or larger 
than the expected statistical error. Therefore, they should be taken into account when 
searching for effects of a strongly interacting scalar sector. In fact, being negative, the 
corrections increase the sensitivity to this kind of effects, which typically appear as an 
enhancement of the WW-scattering cross section. 
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A On-shell projection 

In this section we give the explicit form of the on-shell projection of the two incoming 
W bosons which induce the VBS. When performing the projection, care should be taken 
that the on-shell projected W-boson momenta lie in the physical phase-space region. This 
can be ensured by fixing the angles of the two incoming bosons while performing their 
on-shell limit. We thus fix the direction of the incoming W + in the CM frame of the 
W-boson pair. In this way, the on-shell projected momenta can be written as 




(A.l) 



where \/§ is the VBS CM energy, (3 = \j\ — 4M^/s, and q + is the three-momentum of 
the incoming (off-shell ) W + boson in the CM frame of the W-boson pair. The on-shell 
momenta are then boosted to the laboratory frame. 
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B Logarithmic electroweak corrections 



In this appendix, using the general results of Refs. [12,14], we derive analytical formulas 
for the logarithmic electroweak corrections to the subprocess 

w^fcOwJ, (k 2 ) -> wz X3 (-h)w\(-h). (B.l) 

The momenta fc 3 , k± and the helicities A3, A4 of the final states are defined as incoming. 
The 2^2 process (B.l) is thus equivalent to the 4^0 process 

WM(h)W+ 2 (k 2 )Wl 3 (k 3 )W; 4 (k4) - 0. (B.2) 

This convention facilitates the application of the formalism of Refs. [12, 14], which is based 
on n — > reactions. We consider the limit of high energies and large scattering angles, 
where all invariants are much larger than the electroweak scale, 

\nj\ = \(ki + kjf \ ~ |2£^| > A4 for i ^ j. (B.3) 

In order to specify our conventions for the gauge-boson helicities, which we denote by \ = 
0, ±1, we choose the CM frame. There, the gauge-boson momenta can be parametrized 

as 

fc? = £7(1,0,0,1), -kg = £7(1, sin 7?, 0, cos 7?), 

£# = £7(1,0,0,-1), —k% — £7(1, - sin i?,0, - cost?), (B.4) 

and for the Mandelstam variables we have r 12 = r 34 = 4£7 2 , r i3 = r 24 = — ri 2 (l — cos-#)/2, 
r 23 = ri4 = — 7-12(1 + cos$)/2. Note that mass terms are systematically neglected in the 
high-energy limit. The polarization vectors for transverse gauge bosons (Aj = = ±1) 
read 

e^h, n) = -^=(0, 1, ni, 0), e^i-h, -r 3 ) = -^(0, costf, r 3 i, - sintf), 

e»(k 2 ,T 2 ) = -L(0,-l,r 2 i,0), e^(-k 4 ,-r 4 ) = ^(0, - costf, r 4 i, sin 7?). (B.5) 

Here and in the following, the symbol Tj = ±1 is used to denote the helicity of transversely 
polarized gauge bosons. Longitudinal gauge bosons (A« = 0) have to be related to cor- 
responding would-be Goldstone bosons using the Goldstone Boson Equivalence Theorem 
(GBET) as discussed in Appendix C. There, we also introduce effective couplings for 
longitudinal gauge bosons, which permit to apply the general results of Ref. [12] directly 
to physical matrix elements. 

The matrix element for the process (B.l) or, equivalently, (B.2) is denoted as 

Jvk A i A 2 *3 ^4 = Jvl A i A 2 ^3 *4 (r 12 , r 13 , r 23 ). 

In the high-energy limit, we restrict ourselves to the matrix elements that are not mass- 
suppressed by factors of order M^/^/r^. By means of the GBET, it can easily be seen 
that only the matrix elements involving helicity combinations with an even number of 
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longitudinally (L) and transversely (T) polarized gauge bosons are not suppressed. The 
reason is that all vertices involving an odd number of Goldstone bosons are suppressed by 
coupling factors proportional to masses. In addition, the matrix elements with helicities 
TL — > LT and LT — > TL are suppressed, i.e. 



M w~w+w+w- 



M 



w -w+w+w ( 



0, 



(B.7) 



up to terms of order M^/r 12 . Therefore, in the following we restrict ourselves to the 
non- suppressed combinations 



LL 
TT 
LL - 
TL 
LT 
TT - 



> LL 

> LL 
TT 
TL 

> LT 
TT 



Wo W + • 

w-w+ 

Wo W + • 

w T -w + 
w -w+ 
w-w+ 



w -w + , 
>w -w + , 
w:. 3 w± r4 , 
• w: r3 w + , 
w -wl T4 , 
>w: r3 w± r4 . 



(B.8) 



The amplitudes for LL — > TT, TL — > TL, and LT 
amplitude TT — > LL using the relations 



LT can be obtained from the 



M W ° W W ?3 W TA = M Wr * W ^ W ° W ° 



w-wj,w?w: 



k aW~ W + W + W~ 



j^W W+W+W T4 _ j^w Ti Wr 2 w^ w 



JVl t i T s o vv o 



w-w+w+w,: 



ri2*-*ri3 



(B.9) 



which follow from crossing symmetry. 



B.l Structure of the one- loop logarithmic corrections 

In the following sections, our results for the one-loop logarithmic corrections are given 
either in explicit form, or as correction factors 

wr w+ w+ wr 

5 w- w+ ^w: A3 wl A4 = w _ , (B.10) 

relative to the Born matrix elements. Following Ref. [12], we split the logarithmic correc- 
tions (B.10) according to their origin as 

5 = 8 LSC + 5 SSC + 8 C + 5 PR . (B.ll) 

The double logarithms originating from soft-collinear gauge bosons are split into leading 
contributions <5 LSC of the type a In 2 (|r 12 |/M 2 ) and subleading contributions <5 SSC of the 
type a In (|r 12 |/M 2 ) In (\rij/r 12 \), with = r 13 ,r 2 3. These latter depend on ratios of 
Mandelstam variables, \r 13 /r 12 \ = (1 — cos-#)/2, |r 23 /r 12 | = (1 + cos-#)/2, and thus on the 
scattering angle i? between initial and final states. Purely angular-dependent logarithms 
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of the type a In 2 (1^ /^12 1 ) and a In ( / ^12 1 ) are neglected in our approximation. The 
part 5 C contains the single- logarithmic contributions from collinear (or soft) particles to 
loop diagrams and wave-function renormalization constants. Finally, <5 PR consists of the 
single logarithms that originate from parameter renormalization. 

All these logarithmic contributions depend on various mass scales M=Mw, Mz, M 7 , 
Mh, m t . This mass dependence is separated from the energy dependence by writing 

-mg), (B.i2) 

i.e. we spilt all logarithms into a contribution with scale M-w and a remaining part that 
depends on the ratio M 2 /M^. The logarithms of m t /M^ and Mh/Mw can be found 
in Ref. [14]. Here we give all logarithms of m t /M w as well as the universal logarithms 
of M H /M W , i.e. those that originate from parameter renormalization and collinear sin- 
gularities. Also the subleading logarithms of the type a In (|r 12 |/M 2 v ) In (M^/M^) are 
included. 

In Ref. [12], the virtual electromagnetic corrections have been regularized by an in- 
finitesimal photon mass M 7 = A and split as in (B.12) into a contribution corresponding 
to a heavy photon (A = M\y) and a remaining part which originates from the mass 
gap A <C M w in the gauge sector. The heavy-photon contribution has been combined 
with the weak corrections resulting into the so-called symmetric-electroweak part of the 
corrections. The remaining mass-gap contribution has been isolated into the infrared- 
divergent logarithms L cm (s, A 2 ,m^), /(M{y, A 2 ), l em {m\) which appear in Eqs. (3.7), (3.8), 
(3.10), (3.12), (4.6), (4.7), (4.10) and (4.33) of Ref. [12]. These logarithms contain only 
contributions from virtual photons. 

In Appendix E we provide simple substitutions that permit to generalize the results of 
Ref. [12] to semi-inclusive 2 — » 2 processes, by including the soft-photon bremsstrahlung 
corrections. The resulting logarithms, defined in (E.l), (E.2), and (E.3), are infrared finite 
and depend on the soft-photon cutoff AE. 



Notation 

The coefficients of the various logarithms are expressed in terms of the eigenvalues J^°, 
or of the matrix components /Y",, of the generators 2 

I A = -Q = -\-T\ /* = -£wI + ^7* I^ = Il±lH, (B.13) 
2 c w 2 s w v2s w 

where c\ = 1 — s% = M^/M^- Another group-theoretical object that often appears in 
our results is the electroweak Casimir operator 

r= J2 = + ^T(T + l), (B.14) 

V a =A,Z,W ± c w V 2 / s w 

where T represents the total isospin. 

2 A detailed list of the gauge-group generators and of related quantities that are used in the following 
can be found in App. B of Ref. [12] and App. B of Ref. [14]. 
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B.2 Leading soft— collinear corrections 

Below we list the angular-independent leading soft-collinear (LSC) corrections for various 
polarizations of the gauge bosons. These results are obtained from Eqs. (3.6) and (3.7) 
of Ref. [12] and depend on the eigenvalues of the electroweak Casimir operator 



1 + 2r 2 2 



as well as on the squared Z-boson couplings 



2 / 2 2 \2 

(T z \ 2 — C-w (T z \ 2 — v C w ~ S wJ / R 1R \ 

(v±) - -r, (V) - 4,2 r 2 • ( B - 16 ) 



The explicit expressions for the electromagnetic logarithms L EM (M^ V ), which contain 
contributions from both virtual and real soft photons, are given in (E.l). 
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Mixed polarizations 
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(B.18) 



Purely transverse polarizations 



rLSC 

°w-w+-wz T3 wj 



T 4 



2^ 



C^ln 2 



|Tl2| 



2(4 ± ) 2 ln 



2L* M (M^). 



In 



|TL2| 



(B.19) 



B.3 Sub leading soft— collinear corrections 

The angular-dependent subleading soft-collinear (SSC) corrections to the 2^2 processes 
(B.8) are obtained by applying the formula (3.12) of Ref. [12] [see also (C.2)] to the 
corresponding 4^0 processes which result from reversing the outgoing particles by a 
crossing transformation, i.e. the charges of the outgoing states have to be reversed as 
in (B.6). The corrections to matrix elements involving longitudinal gauge bosons are 



22 



obtained via the GBET applying Eq. (3.12) of Ref. [12] to corresponding matrix elements 
that involve would-be Goldstone bosons. Alternatively, as discussed in Appendix C, one 
can directly apply Eq. (3.12) of Ref. [12] to matrix elements for longitudinal gauge bosons 
using the effective couplings (C.3). 

The SSC corrections originating from soft neutral gauge bosons N = A, Z result in 



E * 



7V,SSC 

w: wi 



N=A,Z Al A2 ~ A3 ~ A4 



_|_ I T N T N -4- T N T N 



■A 4 



In 
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^23 1 
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-4 In 



/-AT T N , j-AT rAf 



■wr 1 w+ 



A2 A4 



In 



ki2| 



rig] 

l r 23| 



/EM 

*SSC; 



;b.2o) 



where the couplings for longitudinal gauge bosons are given in (C.13), whereas for 



transverse gauge bosons = =Fl and = ±c w /s w . The electromagnetic logarithms 
/!£ are defined in (E.2). 



Soft virtual W bosons can be exchanged only in the r 13 channel, and yield 



wr w+ w\ + w; 



v=w± 



X 



a i / r 12 

2~7r [ M? 



In 



w, 



r 12 



. , . , , N Xl W+ N' wr . , . , , Wr" AT A „VF+ AT' 



^2 



;b.2i) 



where the sums over N x and 7V A depend on the polarization A. In the case A = 0, 
they run over N = H,Z , whereas for A = r = ±1 they run over N T = A T , Z T . The 
corresponding couplings , jf and J^., Jf; are given in (C.15) and (C.16), respectively. 
The SU(2)-transformed Born matrix elements that appear on the right-hand side of (B.21) 
are evaluated in the high-energy limit in Appendix D. 

In the following we list the explicit results corresponding to all helicity combinations 
(B.8). Note that, owing to crossing symmetry, both contributions (B.20) and (B.21) are 
invariant with respect to simultaneous exchange of the polarizations Ai <-> A4, A2 <-> A3, 
i.e. 



riV,SSC 

w: w+ ->-w- w H 



E * 

N=A,Z Al A2 



Af,SSC 

wr w+ ->w _ , w+ ' 



E K 

N=A,Z A4 A3 ~ A2 ~ Al 

J2 S v ' ssc M w ^s w ^ = J2 5 v ' ssc M w ^K w m. 
v=w ± v=w ± 



(B.22) 



Therefore, the corrections for LL — > TT and LT — > LT can easily be obtained from those 
for TT -> LL and TL -> TL, respectively. 



Purely longitudinal polarizations 

The contribution (B.20) of soft neutral gauge bosons gives 



E * 



iV,SSC 

wrrw+^w-w+ 



a 

4tt<s^ v c^ v 



In 



In 



r 23 



-4 In 



l r 23| 



1 /EM 

rssc- 



(B.23) 
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The contribution (B.21) of soft W bosons yields 



V 5 V > SSC M W ° W ° W ° W ° = 
v=w ± 



= ^lntel 
2tt lil^ 



In 



F13I 



E J s Jg, 



S,S'=H,Z 



a In SlnM 



\\ -" \\ / \ I ^12 | 



M 



W +Z W - 



Mr° M » + M^ HW ^ H ) - (if - Z ) 



(H <-> Z ) 



(B.24) 



Using our expressions (D.7), (D.8) for the SU(2)-transformed Born matrix elements in 
the high-energy limit, and dividing by the Born matrix element (D.6) we obtain 



ErV.SSC 
5, 



v=w± 



w-w+-w -w+ 



a In W In W 



Ah e 2 r 13 - r 23 
2 2s? 



+ e 



2 2 
2 S w ~~ C w r 12 ~ r 23 



4c2 r 2 



»"13 



A H + 



- r 23 rig - r 23 



'W r 12 
-1 



?~13 



, (B.25) 



where Ah is the scalar self coupling defined in (D.5). 
Mixed polarizations: TT -> LL and LL -> TT 

For the TT — > LL configuration, the contribution (B.20) of soft neutral gauge bosons 
gives 



E & 



iV,SSC 

w-w+^w -w+ 



V=A,Z 

Soft virtual W bosons (B.21) yield 



° In ™ In m 



r 23 



41n O e - (R26) 



2 <^ SSC .M 
y=w± 



r 12|/ at=a,Z S=H,Z 



a .In M In K 



47TSv 



r 12 







s w J 



(B.27) 



Using the Born amplitudes (D.ll), (D.12) in the high-energy limit we obtain the relative 
correction 



\p rV,SSC 

W-W+^W-W + - 2^ Vr23 



° ^-llln 



r 12 



(B.28) 



These results can directly be extended to the LL — > TT configuration using (B.22). 
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Mixed polarizations: TL -> TL and LT -> LT 

For the TL — ► TL configuration, the contribution (B.20) of soft neutral gauge bosons 
gives 



E & 

N=A,Z 



7V,SSC 

w T -w+~»w- W+ 



« i / r 12 

In 



ln(N)+iln /|ri31 



l r 23| 



4c 2 



Soft virtual W bosons (B.21) yield 

y 5 v,^c M w-w+w+w- = ^_ ]n ( hA ln 



r 12 



ri3i 
ki 2 | 



-41n^/ s ™. (B.29) 
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Fl2| 



\ \ri2 / s 



M A T1 W+Z T3 W - + {A ^ 



(H <-> Z ) 



+ 



4 s 2 



M W-HW+H + M W-HW+Z 



(B.30) 



Using the Born amplitudes (D.16), (D.17) in the high-energy limit we obtain the relative 
correction 



E & 



v,ssc 

w^w+-+wz„w+ 



a 



t 3 ,v o 



27TS 2 Vr 23 ' 2 



(B.31) 



These results can directly be extended to the LT — > LT configuration using (B.22). 



Purely transverse polarizations 

The contribution (B.20) of soft neutral gauge bosons gives 



E 8 

N=A,Z 



Af.SSC 

W~W+-<-W" W+ 

r l ^2 — T3 — T4 



7TSf 



In 



ln 



4 m 



l r 23| 



Fl3| 
l r 23| 



/EM 

'ssc- 



;B.32) 



Soft virtual W bosons (B.21) yield 



v=w ± 



-^ln ™ In ^ 

27T \M^J 



E Jn^n' 

ri2 l / N,N'=A,Z 



M n ti w+ n> 3 w~ + ^ jv; 2 7V T4 



= ^ln 



2tt~ V M w 



l r 12| 



^W+A„ff- +M W-A T2 W+A T4 
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(a^*^ + A*^^ 4 ) + {A <- Z) 



„2 



In the high-energy limit, using (D.25) we obtain the relative correction factor 

E 5 w-w+ - W " T w+ T ~ ^72~ ( ) ( O ) ' ( R34 ) 

F=w ± T1 T2 - T 3 -r 4 7TS W ^23 \ iW W/ \l r 12|/ 

B.4 Single logarithms from collinear singularities 

The single-logarithms originating from collinear singularities associated to external trans- 
verse or longitudinal W bosons can be obtained from Eqs. (4.10) and (4.33) in Ref. [12]. 
Here we have included also the logarithms of m t /M w and M H /M W , which can be found 
in Ref. [14]. The resulting corrections to all processes (B.8) can be expressed by the fol- 
lowing general formula, which only depends on the numbers n-r,nL of transversely and 
longitudinally polarized W bosons, respectively. 



<5 W r w + -w-, w+ = E 6 w x . = n T$w T + n i^w , ( B - 35 ) 



i=l 



where 



(B.36) 



the one-loop coefficient of the SU(2) /3-function reads fef^ = 19/(6s^.), and the electro- 
magnetic contributions / EM (M^) are defined in (E.3). 



B.5 Single logarithms from parameter renormalization 
Purely longitudinal polarizations 

The renormalization of the parameters g = e 2 /(4s^c|,) and Ah (D.5) in the Born ampli- 
tude (D.6) gives rise to the relative correction 

tpr ( 5Ah Sg\ Ah m 

Vw ^w -w 0+ = - + [ ^ - ~) (B-37) 

where 

A = + = 2 (j™--l) (B.38) 
ri2 r 13 \r 12 r 13 I 
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in the high-energy limit, and the 't Hooft scale of dimensional regularization has to be 
set to ji 2 = r 12 in the counterterms [12,13]. In the on-shell scheme (including the tadpole 
contributions in the renormalization of Ah) these read [14] 

Sg 

9 



— I -b7v In 



c 2 - s 2 



Ml 



9 + 64-324 
184 



5A H 



At 



Aa(M^) + 




;b.39) 



where b™ z 



,2 Lew I _. 
'w°B + C 



2 Lew 
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384-224)/(644), with the U(l) and SU(2) one- 
loop /5-function coefficients and defined in Ref. [12], and Aa(M 2 v ) represents the 
running of the electromagnetic coupling constant from the scale to M^. Within the 
scheme, these effects of the running are already included in the definition of olq^- Thus, 
the Aa(M 2 v )-terms appearing in (B.39)-(B.41) were not included in our implementation. 



Mixed polarizations 

In this case, the Born matrix element is proportional to the squared SU(2) coupling 

9i 



e /4- A s a result, the parameter renormalization yields 



rPR 

w-w+^w -w+ 



= 5 PR 

w -v 

J 92 



92 



+w: T3 w± T4 
a 
An 



rPR 



w+ u w;rw+- 



9 + 64-324 
184 




+ Aa(M^), 



[B.AO) 



with fefy = 19/(64)- Note that the contribution byft In (\r 12 \/MD from parameter renor- 
malization cancels a corresponding contribution b^ In (|r 12 |/M- 2 v ) in (B.35), which is as- 
sociated to the two transverse gauge bosons (tit = 2). This cancellation is analogous to 
the one observed in Eq. (A. 11) in Ref. [12]. 



Purely transverse polarizations 

Also in this case, the Born matrix element is proportional to the squared SU(2) coupling 
g\ and we have 

5 



rPR 

VV T1 VV T2 —r VV VV 



a 
An 



-b e ^\n 



\ri2\ 



+ 



64 



In 
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+ Aa(M*,). (B.41) 

In this case the contributions In (Ir^l/M^) from parameter renormalization cancel 
only two of the four («x = 4) contributions by$ In (|ri 2 |/M(y) in (B.35). 

C Longitudinal gauge bosons 

In order to apply the general formulas of Ref. [12] to matrix elements involving longi- 
tudinal gauge bosons Vq = W^,Z , these have to be transformed first into correspond- 
ing matrix elements involving would-be Goldstone bosons $ a = tfi^, x by means of the 
Goldstone-Boson Equivalence Theorem (GBET) in its naive lowest-order form 3 , 

J v ^ Pix ...V^...i Pin _ ^1-Q v a J^<pi 1 ...® a ...ip in ^ (C.l) 

where ip^, . . . <p in represent arbitrary particles with arbitrary polarizations and Qya is the 
gauge-boson charge, i.e. Qw± — ±1 an d Qz = 0. 

Particular care must be taken of the angular-dependent subleading soft-collinear cor- 
rections (see Eqs. (3.9)-(3.12) of Ref. [12]) 

= ^EE E A^UQ)'^^^ m 

* n k=ll<kV a =A,Z,W ± \ a, I \l'12|/ I k 

which involve non-abelian couplings lX a lipi that are in general non-diagonal and lead there- 
fore to SU(2)-transformed Born matrix elements with ipy ^ tfi on the right-hand side of 
(C.2). For processes involving longitudinal gauge bosons and Higgs bosons, this formula 
has to be applied to corresponding processes involving would-be Goldstone bosons and 
Higgs bosons. Then the unphysical matrix elements on the left- and right-hand sides of 
(C.2) have to be transformed into physical matrix elements by means of the GBET. In 
doing this, the factors i 1- ^" originating from the GBET (C.l) must be carefully taken 
into account, since the would-be Goldstone bosons and Higgs bosons appearing on the 
left- and right-hand side of (C.2) (and thus the corresponding GBET factors) can be 
different. 

In order to avoid these complications related to the explicit use of the GBET, in the 
following we introduce effective couplings for longitudinal gauge bosons and Higgs bosons, 
which permit to apply (C.2) directly to physical matrix elements. To this end the factors 
ji-Q V a f rom the GBET are combined with the gauge couplings for would-be Goldstone 
bosons and Higgs bosons, resulting into the effective couplings 

TV? - ( _\\l~Q vb jV c -l-Q V a jV- _ jV- -l-Qva jV- _ (_:\1-Q vb tV c 

(C.3) 



3 The relevant quantum corrections to the GBET, which contribute to the collinear single-logarithms 
are already taken into account into the corresponding corrections factors (4.33) in Ref. [12]. 



9 + 64-324, ln 
184 



m t 
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where V c = A, Z, W ± , and $ a 



X are the would-be Goldstone bosons corresponding 



to V a = W^, Z Q . We observe that the relations (C.l) and (C.3) can be regarded as the 
result of a reparametrization of the scalar sector through the the unitary transformation 



•l—Qyu ira 
1 V Q . 



(C.4) 



If one performs this transformation directly at the level of the Lagrangian one can express 
the Feynman rules in terms of the longitudinal gauge-boson fields. This simplifies the 
calculation of matrix elements for longitudinal gauge bosons in the high-energy limit 
(see Appendix D). The Feynman rules for the scalar and gauge interactions of the fields 
Si = Wq, Z q , H and their propagators in the 't Hooft-Feynman gauge read 



ya 



Si! (Pi 



ie/ 5+\ (P2 "Pi) W , 



(C.5) 



s i2 {p2) 




■ 2 f rV a l TV a 2 1 



(C.6) 



s. 
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— i- 



^ 2 ^ 3 S i4 + (2-3) + (2^4) 



(C.7) 



S, 



Hi 



s 



'4 



Sji-P) 



iS 



S+Sj 



f - Ml 



2 ' 



(C.8) 



Here, all fields are incoming, the curly brackets represent an anticommutator and products 
of couplings have to be understood as 



Si Sj 



E 



jV a l jV a 2 



(C.9) 



S k =W±,Z ,H 



The Feynman rules (C.5)— (C.8) are closely analogous to those in Ref. [14]. However here, 
as a result of the unitary transformation (C.4), the hermitian conjugation of the fields 
Vq = Z , Wq generates a minus sign, 



ya- 



—V a 



(CIO) 
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where Vq = < K —iY 1 ~ Q v<^)^+ OI - ; equivalently, W$ = Wq and Z$ = Z . As a consequence, 
the coupling matrices and the Kronecker symbols that involve hermit ian conjugate fields 
in (C.5)-(C.8) have to be understood as 



T s+Sj = ~ I s 1 s^ 6 s+s 3 = -SSiSj (C- 11 ) 
for Si = Wq, Z q . Instead, H + — H — H for the Higgs field, and 



Below we list the explicit expression for all non- vanishing effective couplings (C.3), which 
can be easily derived from Eqs. (B.21)-(B-23) in Ref. [12]. For neutral gauge bosons 
(N = A, Z) we have 

N N A Z 1 — 2s2 

^ww' = &o<y'Iwzi with I W a = —cr, I W<J = <J— , (C.13) 



and 



Iz h = IZz =Snz^— ■ (C.14) 



For charged gauge bosons coupling to S — Z , H we obtain, 



1 



Iwgs — Iswz — ft™' J si Jz — a 7T~i Jr— tt~ ■ (C.15) 



Note that the couplings (C.15) are analogous to the couplings for transverse gauge bosons 
defined in (B.26) of Ref. [12], which can also be written as 

jW"' _ fW- a> _ _r to to _ __ to _ _ C w /p 1R \ 

o 



W 



D Born matrix elements in the high-energy limit 

As input for the evaluation of the angular-dependent subleading corrections (B.21), which 
originate from exchange of soft-collinear W bosons, we need the SU(2)-transformed Born 
matrix elements that appear on the right-hand side of (B.21). The needed amplitudes are 
evaluated in the following in high-energy approximation, i.e. omitting mass-suppressed 
terms. For longitudinal gauge bosons or Higgs bosons, which are denoted as 

Si = W±,Z ,H, S t = W^,Z ,H, (D.l) 

we use the couplings introduced in Appendix C. Transverse gauge bosons are denoted as 

V T a = A T , Z T , W±, V T a = A T , Z T , W?, (D.2) 

with r = ±1. 
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Purely longitudinal polarizations: LL — > LL 

For a generic process 

Si 1 (k 1 )S i2 (k 2 ) -> Sisi-ksjSi^-ki), 
in the high-energy limit we obtain the amplitude 



(D.3) 



V r 12 



+ (2<->3) + (l4->3), 



jy° jV a _ _fi £ £ 4 



where 



Ah 



e 2 M 2 



(D.4) 
(D.5) 



is the scalar self-coupling. 

The Born matrix element for W-boson scattering reads 



M w-w+w+w- = 



e" ( r 13 - r 23 r 12 - r 23 
^H + A o o 1 H 



4s 2 c 2 V r i2 



(D.6) 



The SU(2)-transformed Born matrix elements involving two equal neutral states, S = H 
or Z , give 



M S W ° +SW ° = M™° SW ° +S = e s 



Ah _e___ / rig - r 23 r i3 - r i2 



2 44, V r 12 



(D.7) 



with €h — 1 and e^ = — 1. For the case of different neutral states, S ^ S' with S,S' = H 
or Zq, we have 







-e s e 



1 ^13- r 23 n 3 - ri2 A , si - ci n 2 - r 23 



4s 2 V 7-12 



r 2 3 



+ 



4s?„ 



(D.8) 



Mixed polarizations: TT -> LL and LL -> TT 

For a generic TT — > LL process, 

ow: 2 (^) -> s i3 (-k 3 )s i4 (-k 4 ), 

in the high-energy limit we obtain the amplitude 



M 



= 2e 2 (l -5 n 



T2/ 



r 23 fjV a ijV a z 



r 13 / jV a 2 jV a i 



i 3 *4_ 



Inserting the explicit values of the couplings we obtain 



M 



w-w+w+Wo 



6 T 2 3 

— (1-5 T1T2 ) — , 

S w ^12 



(D.9) 



(D.10) 



(D.11) 
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and the SU(2)-transformed amplitudes 



M N T1 W? 2 HW- = _ M W~N T2 W+H = _ M N T1 W? 2 Z W - 

= — (1 - d T1T2 ) — D N , 



= -M' 



W-N T2 W+Z 



with N T = A T , Z T , and 



ri3 

T23 



1 ri 3 1 - 2sl 



2s w C w T23 2s w C v 



(D.12) 



(D.13) 



Corresponding amplitudes for LL — > TT processes are directly obtained using 



v%v T a 3 3 s i2 s h 



c. c. V S V - 

Mo 12T3T4 =M 



Mixed polarizations: TL -> TL and LT -> LT 

For a generic TL — > TL process 

in the high-energy limit we obtain the amplitude 



-^o =2e 2 (l-5 T1T3 ) 



r 23 / jV a i jV a s\ | r 12 / jV a 3 jV a i 



r 12 / r y«3 r y a 



_Ti3 x /a » 2 fi »4 ^13 - ' "» 2 "*4 

Inserting the explicit values of the couplings we obtain the Born amplitude 



(D.14) 



(D.15) 



si ri 3 



and the SU(2)-transformed amplitudes 



A ^ n 



A ^ n 



2s 2 



,1 5 T1T3 ), 



ri3 



2 s 2 



(D.16) 



(D.17) 



with N T = A T , Z T , /£+ = -1 and /£+ = (1 - 24)/(2s w c w ). 

Corresponding amplitudes for LT — > LT processes are directly obtained using 

c- l/ a 2 e T/ a 4 V a 4 l/ a 2 C 

Vt 2 °«3 Kt 4 _ ^| Vt 4 °*3 Kt 2 °n 

Purely transverse polarizations: TT — > TT 

For a generic process 



(D.18) 
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in the high-energy limit we obtain the amplitude 

U^ 1 ^ 2 ^^ 4 _ o £_ [faSU(2) ,SU(2) _cSU(2) ,SU(2) \ Af'r If \\ 
JV1 * 2 |_\ V a lV a 2 U V a 3V a i U V a lV a 3 V a 2V a 4: J V ' ' I' / 

/cSU(2) rSU(2) _ rSU(2) rSU(2) \r>(->r \\] (y) -\Q\ 

where the matrix # su ( 2 ) is defined in Appendix B.2 of Ref. [14], and has the non- vanishing 

comnonents r^ SU{2) - s 2 ^ SU(2) - r 2 r^ SU(2) - r^ SU(2) - -r s r^ SU(2) - tf SU(2) - 1 
components o AA — s w , o zz — c w , o AZ — o ZA — c w s w , o w+w+ — o w _ w _ — i. 

The functions A, B depend on the invariants and the polarizations f — (ti, r 2 , T3, T4). 

The only combinations of polarizations that yield non- vanishing contributions are 

T a = {a,a,-a,-a), n = (a, -a, a, -a), f c = (-a, a, a, -a) (D.20) 

with a = ±, and we have 

^ fa}) = " — , Mn, {nj}) = A(f c , {r y }) = (D.21) 

?~23 ?"l2?"23 ^12 

and 

B(?a, {r«}) = - — , S(f 6 , {r -}) = £(f c , {r y }) = (D.22) 

ri3 r i2 r 12 r 13 

whereas if f ^ r a , r b , f c then A(t, {r\j}) = B(r, {r^}) = 0. Inserting the explicit values of 
the couplings we obtain 

MF"*" 4 " 5 = {r«}), (D.23) 



and the SU(2)-transformed matrix elements 

A<*>W* = = , {r„», (D.24) 

for iV, iV' = A, Z. We note that the ratio between the Born amplitudes (D.24) and (D.23), 

N^W+K S W T - W~N'W+N Ti 

Mo _ M 1 2 3 _ X SU(2) A(t, { rij }) _ _ A su(2)ri 3 , , 

W-W+W+W- ~ W-W+W+W- ~ N ' N RC-r fr \\ ~ N ' N r ' { > 
^^/^ ^ 1 2 3 4 yV^Q 3 -^ViViji/ r 23 

is independent of the polarizations. 

E Electromagnetic virtual and real contributions 

In this appendix we provide simple substitutions that permit to generalize the results 
of Ref. [12] to semi-inclusive 2^2 processes, by including the soft-photon brems- 
strahlung corrections 4 . These substitutions concern the infrared- divergent logarithms 



4 The soft bremsstrahlung corrections to squared matrix elements factorize into the squared Born 
matrix elements times correction factors. These latter have been divided by 2 and combined with the 
virtual correction factors to (non-squared) matrix elements given in Ref. [12]. 
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L cm (s,A 2 ,m 2 ), /(M^A 2 ), l cm (m 2 k ) that appear in Eqs. (3.7), (3.8), (3.10), (3.12), (4.6), 
(4.7), (4.10) and (4.33) of Ref. [12] and have to be replaced with the logarithms L EM (m 2 k ), 
Zgfc and Z EM (m 2 ), defined in the following. 

These results are valid for arbitrary 2^2 processes in the CM frame, with a soft- 
photon cut-off AE. The contributions from virtual photons (superscript 'em') and real 
bremsstrahlung (superscript 'brems') as well as their sum (superscript 'EM') are given 
separately and split into leading-, subleading-soft-collinear and collinear (or soft) parts. 



Leading soft-collinear contributions 

The terms L cm (s, A 2 , m 2 ,), which are defined in Eq. (3.8) of Ref. [12] and contribute to 
Eq. (3.7) of Ref. [12] have to be substituted by L EM {m 2 k ) = L cm (m 2 k ) + L brems (m|), with 
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(E.l) 



Subleading soft-collinear contributions 



The terms l(M^, A 2 ) in Eqs. (3.10) and (3.12) of Ref. [12] have to be substituted by 



/EM 

f SSC 



/em _i_ /brems ,,.:-(- V, 

'SSC + 'ssc > Wltn 



/em _ <X i ( M W 

/ssc "4^ ln n^ 
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'SSC 



An I AAE 2 ) 
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'SSC 



a , ( M 2 
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An V AAE 2 



(E.2) 



Collinear and soft single logarithms 

The terms Z em (m 2 ), which are defined in Eq. (4.7) of Ref. [12] and contribute to Eqs. (4.6), 
(4.10) and (4.33) of Ref. [12], have to be substituted by l EM (m 2 k ) = l cm (m 2 ) + l hreras (m 2 k ) , 
with 
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(m. 



ibrcms , 
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An 
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kl An 



jEM/ 
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